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Abstract 

The most precise measurements are done at present by timing of ra- 
diopulsars in binary systems with two neutron stars. The timing mea- 
surements of the Taylor-Hulse pulsar B1913+16 gave the most precise 
results on testing of general relativity (GR), finding implicit proof of ex- 
istence of gravitational waves. We show that available results of existing 
measurements, obtained to the year 1993, in combination with the results 
of the Mariner 10 in (1992), give the boundaries for the variation of the 
gravitational constant G/G inside the limits (—0.6 +2) ■ lO"^'^ year~^. 
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Since 1937 [15], [6], and up to present times [18], [7] the possibility of varia- 
tion of fundamental physical constants arise interest among physicists. While 
theories exist where these constants may vary [21], [19], [13], [2], it is evident that 
only from experiments, or better to say from the astrophysical observations, we 
may obtain a realistic answer. The spectral observations of distant quasars have 
been used for estimation of upper limits for the variations of the fine structure 
constant, electron and proton masses [19]. 

Variations of the gravitational constant G have been measured by different 
methods [18], among which are investigations of stellar and planetary orbits. 
The measurements on the Viking Lander ranging data gave the limits [8] for 
the variation of the gravitational constant in the form 

G/G=(2±4).10-iV"\ (1) 
and the combination of Mariner 10 and Mercury and Venus ranging gave [1] 

G/G =(0± 2) -lO^^V^S (2) 

Lunar Laser Ranging experiments are lasting for many years, and improve the 
estimations of G/G from G/G < 6 • 10~"'^^yr~^ in [5] the year 1994, up to 
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G/G= (0.4±0.9) • 10-^V"\ 



(3) 



in [20] the year 2004. 

Recycled pulsars give now the best available timing precision. Especially 
important results may be obtained from timing measurements of the pulsars in 
close binaries consisting of two neutron stars. These measurements had been 
already used for checking general relativity [16] (PSR B1913+16), and [12], [10] 
(PSR J0737-3039). From the long time observations of the the pulsar PSR 
B1913+16 changes in the binary period had been measured, and explaining these 
changes by emission of gravitational waves, it was found [17], that "Einstein's 
theory passes this extraordinary stringent test with a fractional accuracy better 
than 0.4%" . Similar accuracy was obtained from only one year observations of 
the close binary consisting of two pulsars PSR J0737-3039A and B [12] [10]. 
Farther observations of this system should considerably improve the precision, 
in comparison with the one reached by PSR B1913+16. 

The restrictions on the gravitational theory, and checking of its validity is 
possible because timing permits to measure several quantities what gives an 
overlapping information. With the two masses as the only free parameters, 
the measurement of three or more post-keplerian parameters over-constrains 
the system, and thereby provides a test ground for theories of gravity. In a 
theory that describes a binary system correctly, the post-keplerian parameters 
produce theory-dependent lines in a mass-mass diagram that all intersect in a 
single point [3], [10]. Measurements of 5 relativistic parameters together with 
independent determination of the mass ratio, due to presence of the second 
pulsar, are available in PSR J0737-3039, in comparison with only 3 parameters 
in the PSR B1913+16. 

The variation of fundamental constant could be incorporated in the general 
procedure of the timing data development. With 6 measurable quantities at 
present (and possibility to increase this number in future [10]) variability of 
several fundamental constant can be investigated. 

Nevertheless, even without any complicated calculations strong restrictions 
can be obtained to the variation of the gravitational constant using only the 
measurements of the variations of the binary period. Change of G has the most 
evident influence on the orbital motion, because only emission of the gravita- 
tional waves compete here with the variation of G. Influence of G variation of the 
pulsar period changes, star contraction at increasing gravity, or star expansion 
at its decreasing, accompanying by corresponding period changes at constant 
angular momentum, are strongly contaminated with the intrinsic changes of the 
rotational period, due to losses of the rotational energy and angular momen- 
tum. A theoretical description of these losses is rather poor [14] . If we accept 
the correctness of general relativity inside the observed accuracy, than we can 
use the residuals (error box) of the measurements of decay of the binary period 
Pf, for the estimation of variations of G. Such estimations have been done first 
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in [4], and later in [9], using timing data from [16] with the result 



G/G = (4 ± 5) • lO-^Vr'^ (4) 

Here we make slightly more careful analysis of the data in [16], and combine it 
with Mariner 10 data, to decrease the error box for G/G, in comparison with 
[9]. Influence of gravitational waves, and G variations on the binary period 
may be considered separately in linear approximation due to smallness of both 
effects. The orbital angular momentum is conserved during slow G variations, 
as an adiabatic invariant, so only total energy of the system is changing due 
to change of the gravitational constant. In newtonian approximation the total 
energy E and orbital angular momentum M of the binary in the mass center 
system are written as 

(a) E = , (6) M = mav^a- (5) 

Here mi, m,2 are masses, considered as points, m = mi+in2 ' ^'"t"^ ^^'^ azimuthal 
velocity at apogee, where radial velocity is zero. The parameters of the binary 
orbit, such as large semi-axis a, eccentricity e, and orbital period Ph, are uniquely 
determined by the conserved values < and M as [11] 



Gm^m2 , . „ ^ M 2\E\M^ 

(^) OIZ7I ' W A = 7rGmim2W-— -r, {f)e = Jl 



2\E\ ' " "^"^'"^"Y 21£|3' ^ G-^mmlml 

(6) 

Slow variation of G imply the variation of Pp, a, v^ai E and e, which are 
connected by relations following from variations of corresponding relations in 
(5), 6) 

SvAa 5a 5\E\ 5G 5a 

(a) ^li^l = (mvl - ^) ^ + (7) 

m ^ <5G _ 3^ de ^ \E\ fSG _ 6\E\\ 

^ ' Pb a 2\E\' e G^mmiml \ G \E\ J ' ^ ' 

It follows from (7) 

5a ^ _6G Sv^^SG S\E\_^ SG 
a G' v^a G' [i^l G' 

what, using (8), gives 
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Table 1: Error budget for the orbital period derivative, in comparison with the 
general relativistic prediction, from [16]. 

Parameter (10^^^) 
Observed value -2.4225±0.0056 
Galactic contribution -0.0124±0.0064 
Intrinsic orbital period decay P°^^ - if"' -2.4101±0.0085 
General relativistic prediction P^^ -2.4025±0.0001 



5Pt JG G IPb 6e 

^ = -2^, ^^~2n' - = ^'^ = ^-^^^- (10) 

The first relation of (10) was derived in [4] in a more complicated way. So, 
from measurements of changes of the orbital period we may estimate variations 
of the gravitational constant, after account of gravitational wave reaction. For 
the pulsar PSR B1913+16 the error budget for the orbital period derivative, 
in comparison with the general relativistic prediction, is given in Table 1, from 
[16]. 

It follows, that the observed values coincide with GR prediction within the 
error bar 

(^pobs _ pgalypGR ^ 1.0032 ± 0.0035. 

Supposing that all deviations from GR arc connected with G variation, we 
obtain the upper limit for these variation as follows 

g = (4.3±4.9)-10-iVr-^ (H) 

To obtain this result we have assumed correctness of GR prediction (last line in 
the Table 1), used the difference between the observed (third line in Table 1) and 
theoretical values (—0.0076 • 10~^^) for the result, and the combined (theoretical 
+ experimental) error (±0.0086 • 10~^^) for the error box. To transform the 
nondimensional value of P into prediction for ^ in (11), we should multiply 
the result and the error by = 565.4 yr"'^, where Pt = 0.323 days = 8.84 • 
10~^ years. So, more accurate estimations (11) slightly improve the values 
(4), obtained in [9]. This improvement narrows on 40% the boundary in the 
negative region (most important from the theoretical point of view [21], [13], [2]), 
from -1 • 10-12 to -0.6 • IQ-^^. 

Combination of two results (11) and (2) permits to narrow the range of G 
variations, which is now may be situated in the limits 

G/G is within the interval (-0.6, 2) • lO'^^yr'^ 
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These limits are slightly shifted from predictions of superstring theories [21] 
G/G w —1 • 10~^^^^yr~^, and permit for G to remain constant. They are 
independent on the influence of tidal interaction between Earth and Moon. 
Regretfully, the data of measurements of Pb in PSR 1913+16 since 1993, as well 
as the refined data on the PSR J0737-3039 are not available, and the published 
[10] errors in measurement of P;, in the last case arc still considerably larger 
than that for PSR B1913+16 in 1993. Farther observations of both pulsars 
could improve the precision in the estimation of G variations. 
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